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ABSTRACT 

A method has been developed for calculating the internal temperature distribu- 
tion in an aircraft tire while free rolling under load. The method uses an approxi- 
mate stress analysis of each point in the tire as it rolls through the contact 
patch. From this stress change, the mechanical work done on each volume element 
may be obtained and converted into a heat release rate through a knowledge of 
material characteristics. The tire cross-section is then considered as a body with 
internal heat generation, and the diffusion equation is solved numerically with 
appropriate boundary conditions of the wheel and runway surface. Comparison with 
data obtained with buried thermocouples in tires shows good agreement. The 
data presented in this paper were excerpted from reference 1. 
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The origin of aircraft tire heating lies in the hysteretic loss characteristics 
of polymeric materials used for both the nylon textile cord reinforcement and for 
the rubber components in an aircraft tire. Figure 1 shows a typical stress- 
strain curve between the two end points of the cyclic stress states associated with 
rolling a material element through the contact patch of the tire and out again. 
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Figure 1 
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The stress excursion is presumed to occur between two points A and B located 
180° apart on the rolling tire. One extreme is the upper point A where only in- 
flation pressures act to produce an axisymmetric stress state of the tire. The 
other extreme point is the center of the contact patch, where again an axisymmetric 
solution is used to obtain the stress state in the tire deformed against a flat 
surface (fig. 2). 
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Figure 3 illustrates the assumed conditions in the vicinity of the contact 
patch of the tire. A total deflection A is presumed to take place between the 
outer surface of the tire crown and the ground. Inside this contact patch region 
the tire is presumed to be in complete contact with the flat runway surface. 



(a) Assumed geometry at point A 
of figure 2. 


(b) Assumed geometry at point B 
of figure 2. 


Figure 3 
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Figure 4 shows the assumed geometry of the tire at the two extreme condi- 
tions A and B illustrated in figures 2 and 3. In the case of axisymmetric infla- 
tion loads the geometric midline of the tire is represented by a torus as shown. 
For that portion of the tire in the center of the contact patch, the shape is 
considered to be flat in the contact region and made up of a single circular arc 
outside the contact region. The stress solutions are both obtained from axisym- 
metric considerations. 
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The tire is divided into five segments for purposes of both stress and thermal 
analysis (fig. 5). Segment I is that portion of the tire carcass lying in the contact 
patch. It has plane orthotropic elastic properties. Segment II is isotropic 
tread material in contact with the runway surface, and hence is acted upon by 
normal pressures as well as by membrane and bending strains. Segment III is that 
part of the side wall acted on by internal pressure, and also has orthotropic elas- 
tic properties. Segment IV is the sidewall cover of the tire and is considered to 
be isotropic in its material properties. Segment V is the bead region which is 
subject to a special description. 


Edge of 



Figure 5 
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Figure 6 shows the method of splitting the tire cross section into various 
rectangular segments for purposes of thermal analysis. In general the decomposi- 
tion process follows the natural boundary lines of the tire between the orthotropie 
and isotropic regions. 
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Figure 7 is the general heat diffusion equation which is being modeled on a 
discrete basis in the numerical calculation described here. However, the elements 
which are in contact with the outer surface or with the bead area of the tire have 
additional terms associated with either thermal conduction into the metallic wheel 
or convective heat transfer loss to the air. 


K V 2 0 
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dt 
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where 

K = thermal conductivity 
p = density 

o 

q = heat generation rate 
0 = temperature 


Figure 7 
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The equation for convective heat transfer, temperature buildup and internal 
heat generation is discretized into the form shown in figure 8 for numerical 
calculation throughout the tire cross section. Note that due to the higher ro- 
tational speeds of the tire compared to thermal diffusion times, the temperature 
distribution is completely axisymmetric. 


At" “ ^ + a K- 1 , K , l ( 0 K-1 , £” 8 K, 

+ a K,K+l,£ ( 0 K+1, £ -6 K, 

+ a K, 1-1,1 ( 0 K, f,-l -0 K, 9? 

+ a K , Si , £+1 ( 0 K, S,+ 1 _0 K, 

with 

KA 

a i = J— 

l,m,n d 

K = thermal conductivity 

A = area of contact between subscripted elements 

d = distance between centroids of subscripted 
elements 

NOTATION FOR THERMAL ELEMENTS 
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The strain energy density must be related to the hysteretic loss, which results 
in internal heat generation. Due to the lack of appropriate material constants, 
this is now being done on a provisional basis using the loss tangent tan 6, a quan- 
tity usually associated with uniaxial sinusoidal loading. Due to the linear rela- 
tionship between rate of heat generation and frequency, the rate of heat generation 
is linearly proportional to aircraft velocity (fig. 9). 


a ~ n m UV 0 tan6 3 

q = 0.01 r - A/ , 3 cal/cnr/sec 

where 

r 0 = outside radius of tire, in. 
A = tire deflection, in. 

V 0 = aircraft velocity, ft/sec 
U = elastic energy, in.-lb/in.3 

Figure 9 
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In data shown in subsequent figures, calculations are carried out comparing 
predictions of this program with measured temperature rise data. These predictions 
are based on best estimates of material properties taken from the literature. 

(ref. 2). These are given in figure 10. 


MATERIAL CHARACTERISTICS OF 
22 x 5.5 and 40 x 14 AIRCRAFT TIRES 


Region 3 

Young ' s 
modulus 
of rubber , 
E, psi 

Shear 
modulus 
of rubber , 
G, psi 

b tan5 0 , 

psi 

Cord 
angle, 
a, deg 

Thermal 

conductivity, 

K, 

cal-cm/°C-cm -sec 

Density, 
P. deg 

Specific 

heat, 

Cp> 

Tread (II) 


335 

0.15 


5 x 10 -4 

1.0 

0.5 

Carcass 

(I, HI) 

(c) 


0.15 

(d) 

5 x 10‘ 4 

1.0 

0.5 

Sidewall 
rubber (IV) 


335 

0.15 



1.0 

0.5 

Bead (V) 

15 x 10 6 


0.03 






a See fig. 5. 

^Values of tanS Q are given for 25°C. For higher temperatures, tan6 is calculated 
by the approximate expression 

tan6 = tanS o (e"- O1A0) 

where A0 is the temperature rise in °C above ambient. 
c See reference 1. For shear modulus G of rubber use 335 psi. 

^Calculated by cosine law using a = 35° crown angle. 


Figure 10 
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Figure 12 shows comparison between temperature rise and time as measured by 
thermocouples, along with calculations for similar locations on the tire using the 
computer program previously described. Stations A through F correspond to those 
shown in figure 11. 


TIRE:22x5.5-8PR 
SPEED:20 MPH 
FREE ROLLING 
UM TEST (REF. 1) 


O INSIOE 
A MIDLINE 
□ OUTSIDE 
— CALCULATED 


SURFACE: 120 IN. DIA. DRUM 
F z : 4350/4350 = (00% RATED 
8z : 1.341 / 1. 344 = 99.8% RATED 
P 0 : 125/115 =109% RATED 



Figure 12 
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TEMPERATURE RISE,°F 


Figure 13 shows similar temperature vs. time information for a different form 
of this 22 x 5.5 tire. In this case the tire is a 12 PR tire and requires a dif- 
ferent geometric description due to the thicker carcass. This in turn results in 
considerably higher temperatures being generated as this tire is run under rated 
load conditions for approximately 150 seconds. Again, measured data is compared 
with calculations from the computer program. 


TIRE i 22x5.5 12 PR O INSIDE SURFACE M20 IN STEEL DRUM 

SPEED .’20MPH A MIDLINE F z = 7100/7100 100% RATED 

FREE ROLLING D OUTSIDE 8 Z = 1.304/1.312 99% RATED 

UMTEST(REF. 1) CALCULATED Po > 209/190 110% RATED 



TIME, SEC 


Figure 13 
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TEMPERATURE RISE 


Figure 14 shows a cross section of a 40 x 14 28 PR tire used as a model of a 
large aircraft tire. This is a common commercial tire, and in this case a larger 
number of discrete sections is chosen for temperature calculation purposes. Again, 
thermocouple locations are shown by the black dots. Extensive experimental work 
has been done on this tire in order to compare it with the calculation. 


h 

o 


TIRE: 40xl4-22PR 
SPEED: 20 MPH 
FREE ROLLING 
UM TEST (REF. 1) 


O INSIDE 
A MIDLINE 
a OUTSIDE 
— CALCULATED 


SURFACE: 120 IN. DIA. DRUM 
Fz ; 23,500/25,000 = 94% RATED 
8z : 3.378 / 3.378 = 100 % RATED 
P 0 : 160 / 155 • 103 % RATED 



TIME, SEC 


Figure 14 
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Figure 15 shows a comparison between a 22 ply rated version of the tire shown 
in figure 14 and calculations carried out using the program described previously. 
These calculations generally agree will with measured data. 



Dimensions in inches 

Figure 15 
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Figure 16 shows temperature rise data up to approximately 150 seconds of run- 
ning time under essentially rated conditions for the 28 ply rated version of this 
tire. In this case calculations agree well with experimental measurements, although 
here the thicker section tire also requires a separate geometric description. 


TIRE: 40xl4-28PR O INSIDE SURFACE: 120 IN. DIA. DRUM 

SPEED: 20 MPH A MIDLINE F z ; 30,800/ 33,500 = 92 % RATED 

FREE ROLLING ° OUTSIDE 8z : 3.459/3.484 = 99 % RATED 

UM TEST (REF. 1) — CALCULATED P 0 : 207/200= 104 % RATED 



50 100 150 0 50 100 150 c 50 100 150 c 50 100 150 0 50 100 150 0 50 100 150 


TIME, SEC 
Figure 16 
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